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In this letter, we discuss the effect of slow real modes in reaction-diffusion systems close to a 
supercritical Ifopf bifurcation. The spatio-temporal effects of the slow mode cannot be captured by 
traditional descriptions in terms of a single complex Ginzburg-Landau equation (CGLE). We show 
that the slow mode coupling to the CGLE, introduces a novel set of finite-wavelength instabilities 
not present in the CGLE. For spiral waves, these instabilities highly effect the location of regions for 
convective and absolute instability. These new instability boundaries are consistent with transitions 
I to spatio-temporal chaos found by simulation of the corresponding coupled amplitude equations. 
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Amplitude equations have been used in a variety of scientific contexts to describe spatio-temporal modulations of 
reference states close to the onset of criticality. In dynamical systems close to a supercritical Hopf bifurcation [Q, 
the spatio-temporal modulation of the homogeneous stationary state can be described by the complex Ginzburg- 
Landau equation (CGLE). This applies to chemical and biochemical reaction-diffusion systems, among vifhich the 
, Belousov-Zhabotinsky (BZ) reaction is the most well known 
^\ • For chemical systems one of the most striking features of the CGLE is its ability to exhibit spiral wave solutions 
(point defects) similar to experimental observations in a large number of chemical reaction-diffusion systems. Recently, 
it was shown |^ that the CGLE fails to model even qualitatively the dynamics of a realistic 4-species Oregonator 
model ^ of the BZ reaction. This discrepancy is caused by the presence of a slow real mode in the homogeneous 
part of the Oregonator model. By considering a slow-field coupling to the CGLE, we show how the inclusion of an 
amplitude equation for the slow mode gives rise to a finite-wavelength instability for plane waves which is not present 
I ■ in the CGLE. For spiral waves, we have calculated new boundaries for convective and absolute instability. 
O \ Here we consider reaction-diffusion systems whose spatio-temporal dynamics is governed by 



On 
0^ 



ac/9i = F(c;/i) + D-V^c, (1) 



■ where c = c(x, t) depends on the spatial position vector x and time t, and D is a diffusion matrix. Close to the 
k>( \ onset of a supercritical Hopf bifurcation of a homogeneous solution of Eq. (|^) , the spatio-temporal modulation of this 
j_] ■ state can be described by the complex Ginzburg-Landau equation (CGLE). In dimensionless form, the CGLE can be 
written compactly as 

w = w - (1 +iQ;)w|wp -f (1 +i/3)V^w. (2) 

As shown in the two real coefficients, a and /3, and the transformation from w to chemical concentration c can be 
derived rigorously from the original reaction-diffusion system (|l|). The CGLE admits plane wave solutions of the form 
w{t,'x) = Aexp[i(Qx — tot)], with amplitude A = ^/l — and frequency uj determined by the dispersion relation 
io = (3Q^ + (1 ^ Q^)a (where Q = |Q|). The stability of a given plane wave is determined by the growth rate A(fc) of 
perturbations with k || Q 



A(fc) = -(fc2 + 2i(3kQ + A^) ± + a)2 - P + 2il3kQ + A^y. (3) 

In particular, at the Eckhaus border defined by 

D\\=l + afi- 2(1 + a2)gV(l - Q^) - 0, (4) 

a plane wave with given Q, will be unstable to long-wavelength perturbations if -D||(Q) < 0; finally, all plane waves 
are rendered unstable at the Benjamin- Feir-Newell (BFN) instability M where 1 -I- a/? < 0. 
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For simple oscillatory chemical systems, the CGLE shows an almost quantitative agreement with the spatio-temporal 
dynamics of the actual chemical system [^j. However, for more complicated models of the BZ reaction, we have 
previously described Q how the CGLE fails even qualitatively to model characteristic time and length scales of the 
models. This disagreement is caused by the presence of a slow (near-critical) real mode. To incorporate the dynamics 
of the slow real mode into a description valid close to criticality, one may derive an amplitude equation similar to the 
normal form associated with a fold-Hopf bifurcation for homogeneous systems ^ . In dimensionless representation, 
the amplitude equations becomes 

w = w+{l + ij)wz - (1 + ias)w\w\'^ + (1 + i/3)V^i«, (5a) 
ez ^ Xqz + k\w\'^ + eSV'^z, (5b) 

where w and z describe the complex and real amplitudes of the oscillatory and slow real mode respectively. The 
parameter Aq is the reciprocal timescale of the slow real mode and e describes the distance to the Hopf bifurcation 
point. The resonant nonlinear coefficients ag, 7, and k can be derived by application of classical normal form theory 
We shall refer to the system (^) as the distributed slow- Hopf equation (DSHE). The DSHE may be considered as a 
"normal form" or prototype model for oscillatory reaction-diffusion systems with a slow real mode. It describes, for 
example, a realistic 4-species model for the BZ-reaction very well In the following, we shall use Aq = —3.07 x 10"'*, 
7 = —1.56, K = —3.10 X lO"**, and S = 0.67 as calculated for this model, whereas as, f3, and e are regarded as free 
parameters. The value of e in a realistic experiment is of the order 10~^. 

Observe that the DSHE (||) cannot be rescaled to be fully independent of the distance e from the Hopf bifurcation 
point: Except for a rescaling of the amplitude of w, the DSHE converges to the CGLE when e — *■ or Aq ^ —00. In 
either of these limits, the coefficient as is related to the nonlinear coefficient a in the CGLE by 

as = (a + 7K/Ao)/(l + K/Ao). (6) 

In the adiabatic approximation where either the operating point is sufficiently close to the Hopf point (e small) or 
when the real mode Aq becomes large and negative, we expect the dynamics of the DSHE to be fully described within 
the framework of the CGLE. However, the distance from the bifurcation point where the adiabatic approximation 
holds may very well be extremely small and very likely experimentally unrealizable. 
The DSHE admits a family of plane wave solutions of the form 

w{t,x) ^ Aexp[i{Qx- cut)], z{t,x) = Z, (7) 



where A = ■\/(l — Q^)/(l + k/Aq), Z — —A^k/Xq, and the frequency w given by the dispersion relation w = PQ^+aA^ 
with a given by Eq. (||). 

To investigate the stability of the plane waves (^, we consider the growth rate cr(fc) of longitudinal perturbations 
with k II Q. For the DSHE, an analytic evaluation of the spectrum of eigenvalues requires the solution of a cubic 
polynomial with complex coefficients, and is therefore not suitable for analytic evaluation. Instead, we may apply 
second order linear perturbation theory pO| to obtain a series expansion for the growth rates: For the first order 
correction to the Eckhaus criterion (||), we obtain for the DSHE 

where Dy is the Eckhaus criterion (||) for the CGLE. For Q = 0, we observe that the BFN criterion also holds 
for the DSHE. For the CGLE, all plane waves are long-wavelength unstable when 1 -I- a/S < 0. For the DSHE, 
this no longer holds, since a band of plane waves of finite wavenumber still remain stable at the BFN point when 
e > — Ao(Ao -I- k)/{2k{1 + Pj))- This band of plane waves, however, can become unstable to finite-wavelength 
perturbations determined by the condition 

F|| = Rccr'(fc) = 0, |fc|>0. (9) 

For example, for a homogeneous plane wave (Q — 0), expansion of cr(fc) to lowest nontrivial order in e and fourth 
order in k yields 

aik) = -(1 -f aP)e ± 1 [(1 + a^)P^ + ^ ^i^ + f - ~ k\ (10) 

Ao(Ao + K) 
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For the CGLE, the long- wavelength instability will always be the first plane wave instability to occur for the plane 
wave with Q = 0. However for the DSHE this result no longer holds, since the coefficient of order k'^ in Eq. ( [lO| ) 
changes sign when 

(l + a2)/3Ao(Ao + k) ^^^^ 



2{l + a(3-S){a-j)K 



causing a finite-wavelength instability to take place. A similar result can be derived for general Q. A typical band of 
unstable wave numbers are shown in Fig. |l|. We also observe that Eqs. (||) and ( p^ converges to the corresponding 
stability results of the CGLE in either of the adiabatic limits e — > or Aq ^ —oo. The described scenario is illustrated 
in Fig. ||. In (a), e = 10^"', we obtain qualitatively the same behavior as in the CGLE, where the band of Eckhaus 
stable plane waves vanish at the BFN point. In (b), e = 2.0 x 10""*, a finite-wavelength band now persists above the 
BEN point. In addition, a finite-wavelength instability emerges at two points from the Eckhaus curve (open circles 
on figure). As a increases, the corresponding instability curve exhibits a limit point (black circle); above this point, 
all plane waves are unstable. The variation along the finite instability curve of the marginally unstable wavenumber 
of the finite perturbation is shown in (c). 

Similar to the CGLE, the DSHE admits spiral wave solutions (phase defects) in both one and two spatial dimensions. 
In two spatial dimensions, these may be expressed in polar coordinates (r, 9) as w{r, 9) — A(r) exp(i('0(r) — 9 — Lut)) 
and z(r) = Z(r), where A{r) and '0(^) ^re the amplitude and phase of the spiral wave respectively in the complex w 
component whereas Z{r) is the amplitude of the slow z component. These three quantities must satisfy the boundary 
value problem 

A{0) = ipiO) = Z{0) = 0, (12a) 
lim A{r) = v/(l-g.^)/(l + A^/Ao), (12b) 
lim Z{r) = ~A{r)^K/\o (12c) 

r — ^oo 

where Qs = lim,- ''/''('") is the unique wave number selected by the spiral wave. We now discuss the stability 
properties of spiral wave solutions of Eq. (|l^. In order to compare the results with the properties of the CGLE, we 
shall use the parameter a determined implicitly by Eq. (^ as the free parameter, whereas all other parameters in the 
DSHE (|) are kept fixed. 

For large r, spiral waves resemble plane wave solutions of the form Eq. (|7|), and we may therefore expect that the 
spiral wave stability is governed by the corresponding stability for a plane wave with Q — Qg. The transitions where 
spiral waves become either Eckhaus unstable or unstable to finite-wavelength perturbations are therefore D\\{Qs) = 
or F|| (Qs) — respectively. However, as described in [Q, spiral waves emit plane waves with a nonzero group velocity 
Imda/dk, causing perturbations to drift away; the conditions (|^) and(^ are conditions for convective instability, and 
can therefore only be taken as necessary criteria for instability. To determine exponential growth of a perturbation 
u(x, t) even in a steady coordinate frame, we must evaluate the Fourier integral 

1 /"^ 

u{^,t) = — ufc(0)e"=^+'^«*dfe (13) 
27r J-oQ 

for large t in the saddle-point approximation . The crossing to absolute instability is then determined by the two 
conditions cr'(fco) — and ReiT(fco) = 0. 

For four selected values of e, we have determined the variation in the (a, ^-plane of the DSHE instability thresholds 
for Eckhaus, finite- wavelength, and absolute instability as shown in Fig. ^ To solve the associated highly unstable 
boundary value problem, we have used a continuation tool with support for multiple shooting [p^ . The corresponding 
Eckhaus and absolute instability borders for the CGLE are also shown as indicated by the gray-shaded area. Even 
for e small (a), the Eckhaus threshold deviates significantly from the corresponding CGLE curve whereas the absolute 
instability curve almost coincides with the CGLE result. Note that the finite-wavelength instability does not exist 
for this value of e. However, as e is increased (b-d), the finite- wavelength curve completely determines the onset of 
convective instability, and the Eckhaus curve has been omitted from these panels. We observe that both of the DSHE 
instability limits gradually are shifted to the left of the CGLE boundary. Finally, for e = 10~'^ (d), the instability 
limits lie completely outside of the limits predicted by the CGLE. 

As observed for the CGLE ||ri|,Q, the absolute instability (AI) line, is indicative for the onset of persistent tur- 
bulence in the CGLE. Numerical simulations of the DSHE indicated by dots in Fig. ||b and |^ confirms a similar 
observation: spiral waves are convectively unstable below the AI line and absolutely unstable above the AI line, where 
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a transition to sustained turbulence is observed. A representative scenario for the DSHE close to the AI line in Fig. ||d 
is shown in Fig. ^. 

The results derived in the stability analysis presented for the DSHE (H), show that the presence of a slow mode 
in oscillatory chemical reaction-diffusion systems, can give rise to a finite-wavelength instability of plane waves and 
spiral waves, which does not occur in the CGLE. In a real chemical system, this instability occurs at a value of e 
where the amplitude of the oscillations is just below the limit of detection. So even close to the Hopf bifurcation 
point, this instability completely determines the stability of plane waves (which for the CGLE is given solely by the 
Eckhaus criterion Eq. (^)). As shown in both Figs. || and ^, the finite- wavelength instability has profound effects on 
the location of boundaries for convective and absolute stability for spiral waves, and completely alters the classical 
bifurcation diagram known for the CGLE as e is increased. 

For simple model systems of oscillatory chemical reaction-diffusion systems, such as the Brusselator [|l5| and the 
Gray-Scott model ||l6|, the CGLE provides an almost quantitative description of spatio-temporal structures even at 
quite large distances from the bifurcation point; however, models of realistic chemical and biochemical systems, such 
as the BZ-reaction, the horseradish peroxidase system [0, and glycolytic oscillations [ p^Jl9| all posses one or more 
slow modes, and it is therefore unlikely that the CGLE will be applicable for modeling experimental observations on 
such systems. 

We thank Igor Schrciber for valueable assistance with the calculations of the stability boundaries shown in Fig. ||. 
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FIG. 1. Behavior of the real part of the most unstable eigenvalue cr(fc) for a plane wave solution of the DSHE for three 
different parameter values of the parameter e close to criticality (a = 1.85 and /3 = —1). The slow field coupling in Eq. (^) 
causes a finite-wavelength instability to occur before the onset of long-wavelength instabilities. 
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FIG. 2. (a) and (b) show curves in the (|(5|,q) plane describing Eckhaus (EH) and finite-wavelength (FN) instabilities of 
plan wave solutions of the DSHE (j^ for two different values of e and /3 = — 1. (c) describes the variation along the FN curve 
in (b) of the marginally unstable wavenumber of the finite perturbation. 
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FIG. 3. Parameter diagrams showing the variation of the instability boundaries for spiral wave solutions of the DSHE (^) 
for four different values of the parameter e. The figure shows the dominant boundaries for convective instabilities (dashed 
line), Eckhaus (EH) in (a) and finite-wavelength (FN) in (b-d), together with the boundary for absolute instabilities (Al, solid 
line) . The left and right boundaries of the gray-shaded area indicate the Eckhaus and absolute instability curves for the CGLE 
respectively. The Benjamin-Feir-Newell line (BFN) is also shown. For (b) and (d), small circles indicate points where the 
behavior has been confirmed by direct simulation of the DSHE. 
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FIG. 4. Snapshots showing the behavior of the DSHE near the onset of absolute instabihty corresponding to the two 
parameter points indicated by fiUed circles in the bifurcation diagram in Fig. ^. For a = —0.10, a convectively unstable 
transient ends in a frozen spiral state while a = 0.09 gives rise to persistent spatio-temporal chaos. 
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